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Abstract. We study the dynamical properties of a finite dynamical network 
composed of two interacting populations, namely; extrovert (a) and introvert (b). In 
our model, each group is characterized by its size (Na and Nb) and preferred degree (kq 
and Kb ^ Ka). The network dynamics is governed by the competing microscopic rules 
of each population that consist of the creation and destruction of links. Starting from 



o . 

C/3 , an unconnected network, we give a detailed analysis of the mean field approach which 

C/3 ' is compared to Monte Carlo simulation data. The time evolution of the restricted 

degrees (kbb) and (kab) presents three time regimes and a non monotonic behavior well 
^1 captured by our theory. Surprisingly, when the population size are equal Na — Nb, 

(-H ' the ratio of the restricted degree 9o — {kab) / {kbb) appears to be an integer in the 

f*> I asymptotic limits of the three time regimes. For early times (defined by t < ii = Kb) 

the total number of links presents a linear evolution, where the two populations are 
Ps| i indistinguishable and where 6'o = 1. Interestingly, in the intermediate time regime 

^ \ (defined for ti < t < t2 oc Ka and for which Oq = 5), the system reaches a transient 

stationary state, where the number of contacts among introverts remains constant 
rvQ , while the number of connections is increasing linearly in the extrovert population. 

Cn ' Finally, due to the competing dynamics, the network presents a frustrated stationary 

f~^ ■ state characterized by a ratio ^o = 3. 

o 
o 

^H . 
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^ ■ 1. Introduction 



In the last decade, it became clear that the statistical properties of natural or artificial 
networks differ most of the time from the properties of random graphs. Nowadays, it 
is well known that the formation and the dynamics of complex networks are far from 
random but are governed by underlying organizing principles. Some of the most famous 
networks are: the living cells (formed by chemical and chemical reactions), the internet, 
the world wide web, power grids and social networks. See [H El |3l HI [Sj Ej for review 
articles. In addition, it appears that networks are an interesting support for various 
dynamical processes where the topology should play a crucial role. As an example, the 
random walk [3 13 O [101 El [12] or the process of synchronization in complex networks 
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has been some of the favorite playgrounds of the scientific community. One 
should also mention the spread/diffusion of ideas, innovations or viruses [151 CH] that 
are of particular importance. These dynamical problems have been studied on regular 
[m [18], random [191 [201 [21], small world [221 [231 [211 [251 [261 [23 [2H] and scale free 
networks [291 [30]. 

A challenging task, in the characterization of the network topology is the detection of 
sub-units [311 [321 [331 [311 [Ml [361 [SZ], called clusters, modules or communities. These 
are the signature of the hierarchical nature of complex systems [3H1 [391 SO]- Iii the 
case of social networks [H] and in the framework of epidemiology, the spread of an 
infectious disease within a community and its transmission between communities is 
a key to understand and prevent pandemic outbreaks. Interestingly, while for random 
graphs, epidemiological models present a disease free state below a critical infection rate, 
it has been shown, for scale free network, that any spreading rate leads to the infection 
of the whole network [121 [13]. In an effort to build realistic models, epidemiologist 
focused on understanding the robustness of the results to the model assumptions [B] . 
Finding out the simplest microscopic rules of a social network is a difficult task and is 
strongly dependent of the population under consideration. Networks that include social 
communities, social distance and group structures have been proposed [lH [16]. 
Motivated by the sociological interaction amongst extroverts and introverts, we analyze 
the dynamical properties of a finite size network composed of two populations. In 
an epidemiological context, we believe that the daily network of interactions is not 
necessarily scale free. We propose an alternative dynamics of the preferential attachment 
rule, where the populations are characterized by a 'preferred degree', respectively «;„ 
and Kb for extroverts and introverts. In our model, each individual tries to reach the 
'preferred degree' through the possibility to create or destroy bonds. The succession 
of these microscopic events governs the dynamics of the entire network. Evidently, the 
growth of the network is subject to the constraint imposed by the preferred degrees. Such 
constraints are known to lead to single scale networks [17] . Examples are found in actors 
collaborative networks which present a exponential cutoff of the degree distribution. 
We study the time dependence of such a network, starting with no contacts and 
evolving through birth, growth and relaxation. We find three regimes separated by two 
characteristic times ti and ^2 associated to the preferred degree Kb and Ha, respectively. 
In the early time regime, the two populations are indistinguishable while the second 
time regime is understood as a relaxation towards a transient stationary state. In this 
state, the number of links between extroverts is linearly increasing while the number of 
introvert-introvert and introvert-extrovert links are constant, depending on the difference 
of the population sizes. In the last time regime (^2 < t) the system relaxes towards 
the final, stationary state. Not surprisingly, the competing dynamics between the two 
populations induced frustration in the stationary state. This results simultaneously 
of the creation of cross-population links from the extrovert population and of their 
destruction by the introverts. 
The paper is organized as follows: we start, in section 2, with a presentation of the 



Network Evolution Induced by the Dynamical Rules of Two Populations 3 

model, introducing the probability c to create a link as a function of the system 
parameters. In section 3, we give successively a full and coarse-grained description of 
the dynamics of our network through the master equation defined in the complete and 
reduced configuration space. To pursue our analysis further, we derive, for finite system, 
the mean field equations describing the time evolution of the average degree of each 
population. In section 4, starting from an empty network we study the evolution of the 
network through the early, intermediate and late time regimes. For various parameter 
values, the prediction of the mean field theory (in the limit A^ — t- oo) are compare to 
Monte Carlo simulation data. The limitations of our approach are presented in section 5, 
especially in the intermediate time regime where the system presents unexpected finite 
size effects. We summarize our results in the last section and discuss possible extension 
of our study. 

2. The model 

Generally, a realistic network is composed of an arbitrary number of populations or 
groups. In our work, motivated by the competing aspect of different populations, we 
focus our attention on the interaction between just two groups: extrovert and intro- 
vert individuals. The model under consideration is composed of N nodes divided in 
two populations, respectively a (extrovert) and b (introvert), of size A''^ and A''^ with 
N = Na + Nb being constant. Each group is characterized by its preferred degree «:„ 
and Kb. Extrovert and introvert individuals prefer to have respectively many and few 
contacts, which is modeled by imposing «;„ ^ Hb. It is the succession of local events, 
such as creation or destruction of edges by each individual, that governs the dynamics 
of the entire network. Starting with a totally disconnected network (i.e., zero contacts), 
we study the growth of connections induced by the microscopic rules. Of course, more 
sophisticated versions, including the birth and death process of nodes or "mutation" 
from one group to the other, can be considered. 

In our simulations, one time-step is defined by the successive random selection of 
N individuals where each node is updated once, according to the rules detailed below. 
Defining 

A = {Na- Nb)/N, (1) 

the probability to select a node of the population a is given by (1 + A)/2 (and (1 — A)/2 
for b) such that on average, all nodes are updated once in each time-step. We will label 
a node by j, denote its social group by g(j) (g(j) = a,b), and its actual degree by kj. 
When a node is chosen, we will create or destroy a link attached to it, with respectively 
probability c or 1 — c. For simplicity, we will choose the same structure of the probability 
of creation c for the two groups. The only difference being the preferred degree, Kg(j), 
which enters through the dependence on kj in terms of the ratio 

rj = kj/n^ij) . (2) 
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Clearly, we must insist that c = 1 if the degree kj = 0, since there can be no more 
links to destroy. On the other hand, a natural condition is c — > if kj ^ '^g(j)- The 
crossover from 1 to can be characterized by some width, which we also introduce as 
a parameter: a. We may regard this as a measure of the 'flexibility' or 'tolerance' level 
of the population, with larger a corresponding to higher 'rigidity' or lower ability to 
tolerate a degree k far from the preferred value. Specifically, we choose 

<-) = Y^^ (3) 

and denote, for the node j, 

Cj = c{rj). (4) 

as the probability to create a new link with a randomly chosen individual who is not 
already in contact with j. On the other hand, with probability dj = 1 — Cj, we randomly 
choose one of its existing links to destroy. After the update, its degree k'^ is simply 
kj ±1, a restriction which can be easily lifted in future studies. 

Clearly, form equation ([3]), the 'width' of the crossover region decreases with 1/cr, 
so that c{r) approaches the Heaviside function H{1 — r) in the limit a — )■ oo. Much of 
our analysis will be based on this limit and its associated predictions for the system's 
behavior are expected to be universal. Note that the probability c is the familiar Fermi- 
Dirac distribution, by identifying k,K,a with e^,e'',/3. In our context, it has another 
advantage, namely, the ratio of destruction to creation probabilities is a simple power 
of r: d/c = r'^ . Other choices for c are of course possible. However, we believe the 
detailed specifics of c will lead to only non-essential properties displayed by the system. 
For most of our simulations we choose: A^ = 1000, A = 0, Ka = 100, k^ = 10, and 
various a's in the interval [2, 100]. 

3. Master equations and a mean- field approach 

3.1. Full description of the system 

Such interwoven systems are particularly well described by the powerful approach and 
tools of statistical mechanics. It is clear that only few exact results can be extracted and 
that an analytical description remains challenging. The size of the configuration space 
illustrates perfectly the complexity of networks related problems. We start by defining 
the contact variables f^iij between the nodes i and j, such that fiij = 1 if a link connects 
these points and fiij = otherwise. As we exclude 'self-links' in our model, /ij^j = at 
all times. Thus, a configuration, X, is uniquely labeled by the symmetric (off diagonal) 
'matrix' /ij j. Since there are C = N{N — l)/2 elements in such matrices, there are 2^ 
possible configurations (i.e., X = 1,2, ...2^). The master equation of the probability 
PtiX) to find the system in a particular state at a time t, has the usual structure 

^m(^) = E ^xPtiY), (5) 
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where W^ is the probabihty to make a transition from Y to X. Note that here, t + 1 
is not a time-step as defined above, but a single attempt to update the network. Since 
one hnk is always created or destroyed, the probability for any X to remain unchanged 
is zero. Of course, it is possible to write an explicit expression of the probabilities Wx 
in terms of the elements of the set {cj}. However, this approach is intractable and so, 
we will not pursue it further. We will consider simpler approximations below. 

3.2. Reduced phase space 

Instead of the detailed picture above, we believe most of the essentials of this network will 
be captured by a coarse-grained view, namely, a 3 dimensional phase space spanned by 
a (/9), the total number of links between extroverts (introverts), and x? the total number 
of links between the two groups (see Fig. [1]). Defining Pt{a,x,P) as the probability to 
find the system in a macrostate defined by the triplet [a, x, /?), its master equation takes 
the following form 

Pt+i{a, X, /3) = J2 ^*^" + '' ^' /5)^".=F(a + ^' X) (6) 

j=±i 



j=±i 

Y, Pt{a, X, /3 + i)W^,Tix, /3 + 0, 



j=±i 
where 1^^,+ and VF^_ are the probabilities to create and destroy a link of type u 
(z/ = a,x,P)- Even in this simplified picture, the exact transition probabilities {W^_±) 
during one time step are not easy to write, since that would require the knowledge of how 
the a, X and /3 links are distributed in the system. To proceed, we assume that all the 
nodes of the same population have the same degree, so that only their averages, denoted 
by kg (g = a,b), enter into our W. Clearly, the degrees ka and kb receive contributions 
from links within the groups a and /3 and those across groups (x). Specifically, defining 

^ (7) 

= h'> + J^- (8) 

kg/ Kg and construct the creation probabilities: 

(9) 

Let us emphasize, instead of having A^ individual Cj's, we deal with only two "group" 
c's here: Ca and Cb- Next, noting that the creation of a link by an individual requires 
the random selection of another node not already connected, we will approximate that 
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probability by 

P. = ^^^^^^, S = a,b. (10) 

Thus, within the group a and b, the probabihties Wa^+ and Wi^^^ are given by 

Wa,+ i<^^X) = J^PaCa, and l^/3,+ (X,/3) = ^PfeCfe. (11) 

Note that the probabihties {W) depend on only two of the three variables {a,x,P)- 
Following the same lines and defining dg = 1 — Cg, we arrive at 

W^«,-(«,X) = §^4, and W^/5,-(x,/3) = ^^4, (12) 

N ka N kb 

where k^^/k^ is the probability to pick, among the existing links, a link connecting a 

node of the same population. Finally, for the creation and destruction of cross-links, we 

have 

PV^,+(a,x,/3) = ^ j;^ ^" + A^ N ^^' ^ ^ 

H^x,-(«,X,/3) = ^^4 + ^^4. (14) 

As usual, if Pf(a,x, /3) is known, it can be used to find the averages of all observable 
quantities. A relevant example is given by the average degree associated with each 
population at any time t: 

{k^),= Y.k^P,ia,x,(3). (15) 

For convenience, we will drop the subscript in {•)t below, except when its presence 
is needed for clarity. We will also discuss time independent (or almost so) quantities, 
which will be denoted by a suitable superscript (e.g., (•)*). 

3.3. A mean-field approach 

Even in the reduced phase space, the master equation ([6]) cannot be solved exactly. 
To gain some insight into our system, we proceed with a mean-field approximation 
(i.e., replacing averages of functions by functions of averages). In this process, we will 
often encounter creation probabilities that depend on {kg^ . Therefore, let us define for 
convenience 

f = ,, , _ , — , and h = ,, , _ , — . (16) 

From the exact equation dt{i') = {Wu,+) — (W^-), we find 

rj / \ Na Na — 1 — [kaa) „ Na {kaa) f^ „^ /^ „n 

dM) = ^ ^ / - ^^^(1 - /)> (17) 

which, recalling kaa = '^a/Na and rescaling t — )■ Nt, leads to 

dt{kaa) = (l + A - 2^^^) / - 2^(1 - /), (18) 



Network Evolution Induced by the Dynamical Rules of Two Populations 




Figure 1. Reduced phase space of the network 



and similarly 



dt{hb) 



1 - A-2- 



(W^^_2Mn 



h). 



(19) 



N )~ ~ {h) 

The equation for the restricted degree {kab) is also straightforward, though slightly 
different: 

dt{kab) 

With the equations 



l^_M)„,,)_W,i_;)_M(l_,). po) 



{ka) = {kaa) + {kab) / {I + A), {h) = {hb) + {kab) / (l " A), (21) 

we have a closed set of three ordinary, nonlinear differential equations for three functions: 
{kaa), {kab), {kbb) ■ Of course, we can make linear combinations of the above to write 
equations for {ka) and (fcfe). However, this two equations do not close since both involve 
{kah)- One instructive combination is given by p 



P = -^—{kaa) + -{kab) + -^—{hb) 

■{ka) H — {h) 



(22) 



which is the related to the total number of links in the entire system by the relation 
Np = {a) + {x) + (/3). Its evolution is governed by 

dtp= [{1 + A)f+{1-A)h-1] 



— [(1 + A) (1 + {ka)) / + (1 - A) (1 + {kb)) h] . 



(23) 
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A convenient simplification occurs if we take tlie limit iV — ;■ oo, namely, 

d,{Ka)= (l + A)/-2i^(l-/), (24) 

dtih,) = (l-A)h-2^-^{l-h), (25) 



dtikab) =^-^{f + h)-{kat) 



1-f , 1-h 



(26) 



{ka) {kb) 

In tlie following, we will investigate the time evolution of these quantities, using both 
Monte Carlo simulations and the equations above. 

4. Time evolution 

In this section, we present the studies of how an initially empty network evolves to the 
steady state. Such a scenario may be realized when the population is formed from a 
collection of 'complete strangers,' such as at the start of a cruise or the beginning of 
a school year (for, e.g., kindergarteners or college Freshmen). Simulating the network 
according to the rules above, we discover three time regimes, as well as non monotonic 
behavior in the time evolution of the restricted degree {kab)t and {kbb)t- See figure (|2]). 
As we will show, the different regimes can be understood in terms of the mean-field 
approach. In particular, provided the k's are well separated, the characteristic times 
(ti and ^2) separating these regimes become quite well-defined for systems with large 
a. Such cases correspond to a highly 'rigid' population, for which both creation and 
destruction of links are highly probable even for small values of \{kg)t — ^g|- We will 
analyze the three regimes successively, for a range of values of the system parameters. 

4-.1. Early time regime 

In the early time regime, the network is nearly empty, so that the probabilities of creation 
of links are close to unity. Then, the average number of links (u) [v = a, x, P) grows 
linearly with t and it is easy to find the slopes for each of the three classes of links. In 
the limit a — ?■ 00, the probabilities of creation Cg become Heaviside functions and we 
can define the early time regime by the times t < ti, in which Cg = 1. Equation fl23l) 
simply leads to p{t) = t while fl2^ . fl25|) and fl26|) predict trivially 

(A;aa)t = (l + A)t, {kbb)t = il-^)t, {kab)t = {l-^^)t. (27) 

We can include finite population size effects and use (11811201) to obtain rounding of this 
linear growth. In this regime, we should note that 6'a defined by the ratio {kab)/{kbb), 
is time-independent and given by 

^A = 1 + A, Vt < ti, (28) 

with ^0=1- Remarkably, even if the size of the populations are different (A 7^ 0), the 
growth of the degrees of the two groups are indistinguishable: 

{ka)t = {kb)t = 2t, (29) 
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Figure 2. Comparison of the mean field predictions (lines) with the simulation results 
(symbols). The time evolution of the degrees (fcfcb)(i) and {kab){t) arc plotted for 
different values of a, with Kq = 100, ki, = 10 and A = 0. 



a result that breaks down for t ^ ti and a = O (1). Now, we see the meaning of ti: It 
marks the time when the introverts reach their 'tolerance' level, i.e., {kb)t^ = Kb. Thus, 
we have 

ti = fi:b/2. (30) 

If we include finite A^ effects, then this becomes — "f-ln (l — ]^^)- In figure ([3]) and 
(jl]) these results are compared to the Monte Carlo simulation data. The short time 
evolution of the restricted degrees, (/cgg')t, are plotted for a = 10. The linear increase 
of the degrees is in perfect agreement with our predictions. The comparison for various 
values of A are shown in Fig (j4]). In the latter figure, we plot the time derivative of 
the degrees at time t = {dt{kaa)\t=o, dt{kf,b)\t=o and dt{kab)\t=o) as a function of A. 
Agreement between theory and simulations is again excellent. 



4-2. Intermediate time regime 

Similar to considerations above, let us define ^2 by the condition {ka)t2 = i^a which 
satisfies the extrovert population. In the regime, ti <t <t2, with k^ > Kb (and relatively 
small A), the behavior of the system is characterized by two com]9etm(7 tendencies. The 
extroverts have yet to reach their preferred connectivity and so, continue to attempt 
creating links. Since some of these will be cross-links (x), the introverts find themselves 
having more connections than the preferred level and tend to destroy links. In the case 
f^b ^ i^a , we can expect ti <^ ^2 and the nature of this intermediate 'state' is more 
exposed. In particular, a very large system (e.g., A^ ^ Ha) may reach a 'steady' state, in 
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Figure 3. Time evolution of the average restricted degrees (kaa) (circles), (kbb) 
(squares) and {kab) (diamonds) in the early time regime. The data are extract from 
simulation (symbols) realized on a population of size N ~ 1000, after 10 realizations, 
for Ka = 100, Kb = 10, A = 1/2 and <t = 10. The lines are given by the mean field 
results, equations ((27|). 
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Figure 4. Slope of the average degrees {kaa), (kbb) and {k^) at time t = 0, as a 
function of the parameter A. The data are extract from simulation (symbols) realized 
on a population of size N = 1000, after 100 realizations, for Ka ~ 100, Kb = 10 and 
a = 10. The lines are given by the mean field results, equations (|27| . 



which the extroverts continue to create hnks at a constant rate. On average, a constant 
fraction of these new hnks will be cross-links, which are then destroyed by the introverts. 
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In such a state, the rate of creation of cross-hnks by the extroverts balances the rate 
of destruction by the introverts. Thus, both restricted degree (kab) and {kf,b) remain 
essentially constant, values which will be denoted by (kab)* and (kbb)*- The continued 
growth in (fc^) must be entirely through that in (kaa) (which is ~ (1 + A)t). Finally, 
the end of this regime is marked by the extroverts reaching their preferred Ka- We can 
estimate ^2 through k^ = {ka)t2 with 

«^a = {kaa)t, + |^. (31) 

We will show that (kab)* = O (/€(,), and so, for Kb "^ i^a , we find 

h = Y^ ■ (32) 

In the remainder of this subsection, we will demonstrate that this scenario is borne out, 
at least qualitatively. 

Remaining with the simple (but extreme) case oi a, N ^ oo and approximating 
the creation probability for the extroverts by unity (i.e., / = 1) leads to a consistent 
solution. By ignoring terms with (1 — /), {kaa) formally decouples and we need to solve 
a system of only two variables, i.e., (kbb) and (kab)- Specifically, equations ( IMl) . (l25ll 
and fl26|) reduce to 

dt{kaa) = (l + ^) , (33) 

and 

dt{kbb)= il~A)h-2^il-h), (34) 

{h) 

dtiKb) = ^^(1 + h) + M(i _ h). (35) 

2 {kb) 

Since the dynamics of the extroverts decouples (except through an external, "source- 
like" term in the second equation), it is sensible to consider an equation for (kb) in lieu 
of ([35]): 

..(..)^i±^.^.-(M.l)(l-.) (36) 

-^ + ^'^-\-#^(l-/.). (37) 



2 2 ih) 

Since h is a. function of {kb) alone, the latter equation together with equation ( 134|) form 
an explicitly closed set. At large times, we have the solution 

{kaa) = {l+A)t, {k,,)^{kbb)\ {kb)^{kb)\ h^h\ (38) 

where the superscript (*) denotes (quasi-) stationary values of the relevant quantities. 
In particular, eliminating the last terms from equations ( 134|) and ( l37|) . we easily find 

h^ = i^, (39) 
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and so, for cr ^ 1, 

{hV = ^b. (40) 

Notice that these results are 'universal' in that they do not depend on the explicit form 
of the probability of creation ([3]), but only on its being well approximated by a Heaviside 
function. The interpretation of the latter is particularly appealing. The introverts are 
satisfied, as far as their overall degrees are concerned. The only dynamics left in this 
regime is a redistribution, from cross-links {{kab}) to intra-group links {{khh)). Before 
the (quasi-) stationary state is reached, the increase of cross-links are, of course, driven 
by the extroverts (trackable to the first term in equation (135|) ). 
The other (quasi-) stationary values follow readily: 

{hbV = 1^^^. {k^^r = (5 - A)(W*, (41) 

with the ration 6^ = (kab)* / (kbb)* = 5 — A and ^o = 5. In the lowest panels of 
figure ([2]), for which A = 0, cr = 100 and N = 1000, this (quasi-) stationary state is 
well represented, in the range 50 < t < 80. The solid lines are results of numerical 
integration of the full set of equations (Elj), (l25l) and ( 1261) . They confirm this scenario 
quite well. An example is the ratio ^o = (kab)* / (kbb)* = 5. The dots are results from 
simulation. Clearly, there is also qualitatively good agreement, especially considering 
what a drastic approximation went into formulating this simple theory. 

4-3. Late time regime 

The onset of the last time regime is characterized by the following. While finally the 
extroverts reach (on average) their preferred degree k^, the system is not completely 
"satisfied" as the introverts have been driven, up to this time, by the zealous creation of 
cross-links from the extroverts. Thus, readjustments will take place, through mostly the 
destruction of cross-links. The decrease of the number of cross-links x iii turn affects the 
introverts, who compensates by creating more links. Again, the last panels of figure ([2]) 
showcase this scenario quite well, starting at t ~ 80. With both degrees being mostly 
satisfied, the two populations spend this final time period readjusting the distributions 
of cross-links and intra-group links, until the true stationary state is reached. 

Let us denote the relevant quantities of the final state by a superscript (*), and 
seek the solution to dt{k^^i)* = 0. Again, for simplicity, we present the results in the 
limiting case a,N ^ oo. First, considering ( 123|) . we find 

(l + A)/* + (l- A)/i* = l. (42) 

This equation expresses a universal, and intuitively understandable, aspect of the 
stationary state. Dividing both sides by 2, the left can be recognized as the probability 
to create a link in the system, regardless of which group the individual is a member. 
Therefore, equality ( H2|) simply states that this probability must be 1/2. The other 
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equations for /* and h* can be obtained by rewriting the right hand sides of dt{k^^* = 0: 

(kab)* 



(i + A)r-2 

{1-A)h* -2 



;i+A)(u* 

(kab)* 



(i-r) = o, 

'l-h*) =0. 



(1 - A) ih) 

Since /* and (ka)* are uniquely related to each other (and similarly for h* and {kh) 
one has 



(43) 
(44) 



[i + A){k,y 



(i + A)r 

2(1-/*) 



1 



{ka,r = (1 - A) {kbT 



(1 - A) h* 
2(1 -h*) 



(45) 



which provides the necessary second equation between /* and h*. Note that, like 
equation (H2l) this does not depend on the explicit functional forms of the probability Cg. 
To continue further, we can recall the explicit forms, such as {ka}* = Kq (1//* — 1) ? 
or consider the limit cr —t- oo, in which we simply have, regardless of the details of Cg] 

{k^y = Kg. (46) 

Then, the only parameters which control the fixed point are 

R = Kt/Ka and A = {Na-m)/N, (47) 

which appear as coefficients 

A = (l + A)(3 + A)-i2(l-A)(3-A), (48) 

fi = (l + A)(4 + A)-2-i?(l-A)(4-A), (49) 

C = 2A-i?(l-A), (50) 

in a quadratic equation for /*: 

A{f*f -Bf* + C = 0. (51) 

A similar equation is easily obtained for h*, by replacing (-R, A) — )■ {R^^,—A). The 
solutions are readily found, though the expressions do not convey clear messages. To 
gain some insight, we consider two simple limits: R = 1 and A = 0. 

- Of course, the physics of i? = 1 is trivial, since the system consists of two identical 
populations (k^ = i^b = ^)- Thus, (fcg)* = k and f* = h* = 1/2. Here, A represents an 
arbitrary partition (perhaps based on other attributes) of a homogeneous population. 
We find easily (kab)* from considering the fraction of the average degree (ka)*, say, that 
are connected to the "introverts," i.e., {ka)*Nb/N = {ka)*{l — A)/2. But, from equation 
(IH]) this fraction is {kab)* /{I + A), so that 

{kab)* = <l - A2)/2. (52) 

Naturally, this is also the result if we consider the fraction of {kh)* connected to the 
extroverts. From this, it follows immediately {kaa)* = {kbb)* = ^^(1 + A^)/2. 
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- On the other hand, if the two populations are of the same size (A = 0) but distinct 
(na 7^ Kb), the results are slightly less trivial: 

{kaa)* = t^a- {kab)* , {hb)* = Kb - {kab)* , (53) 

1 



{kabY 2 



Ka + Kb \/ ^a l^aKb + Kj^ 



(54) 



An interesting, and somewhat unexpected result is obtained for the limit Kb ^ Ka 
{R <^ 1), in which the average number of links within the group of introverts, (kbb)*, 
approaches a finite fraction of Kb. From the above, we have 

{kabr = lKb{l + 0{R)), (55) 

so that only a quarter of an introvert's connections are with other introverts. Again, the 
ratio (kab)* / (kbb)* appears to be an integer: 6'o = 3. The dynamics in this scenario may 
be referred to as "frustrated" and described as follows. In their attempt to establish 
large number of contacts, the extroverts are just as likely to hit upon the introverts and 
so, to increase kab- But this activity is not welcome by the introverts and links would be 
cut in response. If an introvert were connected to just as many cross-links as intra-group 
links, then both kab and kbb are equally likely to decrease. But such a situation cannot 
be sustained, since only the cross- links are continually being formed (by the extroverts). 
Thus, there must be more cross-links for a sustainable balance to emerge. Examining 
the last panels of figure ([2]) again, we see that this ratio is quite consistent with the large 
t behavior of simulation data. It would be useful if a simple and intuitive argument can 
be advanced to arrive at the ratio 6*0 = (kab)* / {kbb}* = 3. 

With a good grasp of the final stationary state we obtain a general picture of the 
relaxation dynamics through the standard linear stability analysis. Defining the vector 
(j)(t) = {6aa, ^ab, ^bb), with 6gg' = (Ajgg/) — {kgg')* , oue easily obtain the expression of the 
matrix M that governs its evolution: dt4>{t) = —M.(f){t). Three relaxation times emerge, 
given by the inverse of the eigenvalues of M. In the limit a >> 1, a careful examination 
reveals that the largest relaxation time r, which governs the asymptotic behavior, is of 
order one while the two remaining characteristic times scales as 1/a. The details of this 
analysis can be found in appendix. With more effort, we can compute the precise value 
of r and compare the e~*/^ tail to those in simulations. However, we doubt that such 
an exercise will provide any further insight, since we expect the agreement between a 
mean field approach and the fully stochastic problem to be good only at a qualitative 
level. 

5. Limitations of mean field theory 

From the previous section, we see that the simple minded, mean-field approach (with 
large a, A^) provides a good qualitative picture of the evolution of our network. On 
the other hand, it is also clear that there are limitations, which are not easily 
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anticipated. This section is devoted to a brief description of some of the shortcomings 
and speculations on the missing crucial ingredients. 

5.1. Drifts in the intermediate regime 

From our simulation data, with cr ^ 1 and A^ = 1000, ..., 5000, we observe a slow drift 
- that scales with 1/A^ - in all quantities associated with the introvert population, e.g., 
(kab), (kbb), etc. Most significantly, both (kb) and (cb) drift with t. But, in a mean 
field approximation, the latter probability (which is ( [1 + {kb/nb)"] ) ) is replaced 

hj h = [1 + {{kb) /nbY]' and so, becomes intimately related to {kb). In particular, 
if 0" ^ 1, small changes in {kb) would induce drastically magnified (or suppressed) 
changes in h. Examining figures |5] and El it is clear that the drifts in these quantities 
are comparable, rather than dt{cb) ~ adt {kb) (or dt{cb) = 0). Such results are clear 
signals that correlations like {cb (kb)) — Cb {{kb)) are non-trivial. Even if we ignore such 
difficulties and examine the O {1/N) corrections in the mean field equations (ll8H20p . we 
can show that they cannot produce a consistent set of results that includes a linear t 
dependence with O {1/N) coefficients. A more interesting behavior is the scaling of the 
restricted degree {kbb) with the size N of the population. By choosing the origin to be 
ti and using the scaled variable {t — ti) /N, the data for all the A^'s can be collapsed 
(shown in the inset in figures [5] and [6]) . Given that mean field theory cannot even predict 
a linear drift, we believe that an explanation of scaling can only emerge with a better 
understanding of the stochastics. To initiate further studies along these lines, one of us 
(RZ) has investigated the distribution of degrees for the simplest version of this model 
possible - single population, c{k) = Q {k — k) - and found it to be an exponential [18] . 
We believe that such non-trivial distributions will provide some the crucial ingredients 
needed to establish a theory for both the existence of these slow drifts and their scaling 
functions. 

5.2. Finite a effects and universality for large a 

With moderate a's, the creation/destruction probabilities are rounded, so that mean 
field theory predicts smoother transitions between the regimes. This rounding is 
especially pronounced for the second transition, t ~ ^2 • As can be seen in figure ([2]), this 
transition is obscured even for a = 10, displaying only smooth, non-monotonic behavior. 
In other words, finite o rounding in the theory becomes significant when a drops below 
50. By contrast, the noise in the simulations completely masks this transition, even for 
a = 100. At the same time, (apart from a constant off-set) the three sets of data with 
cr > 5 are essentially identical through this transition. Again, we view such differences 
between predictions and simulations as limitations of the mean field theory. 
Another interesting aspect of finite a is 'frustration.' Since everyone is more 'fiexible,' 
the extroverts/introverts can tolerate a finite level of cross-links being destroyed/created 
by the other group. As a result, (r^)* = {ka)*/Ka < 1 and (rf,)* = {kb)*/Kb > 1, as seen 
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in simulations. Naturally, A = {Na — Nh)/N also plays a role. In general, the theory 
above predicts /*, h* ^ 1/2. Given (r^)* = (1//* - 1)^^'^ etc., we see that the ratios will 
differ from unity for finite a. We can gain some insight by rewriting equation (H2l) in a 
manifestly symmetric form: 



[((ra)T-A][((r,)T + A] = l-A^ 



(56) 



Clearly, the two ratios lie on opposite sides of unity. By computing their precise values, 
we see that they agree reasonably well with data (e.g., for t = 300 in figure [2]). In this 
connection, we note that the agreement is much better for intra-group links (e.g., (kbb)*) 
than for the cross-links. Presumably, this aspect is also a limitation of mean field theory. 
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Figure 5. Average degree (fc^f,) as a function of t for different system sizes A^ = 1000 
(circles) , 2000 (squares) , 3000 (diamonds) and 4000 (triangles) . The data are extracted 
from simulations realized on 500 realizations, for A = 0, Kq = 500, Kb = 10 and 
a = 100. Inset shows data collapse of (fcbh) and (fcafc) when plotted against [t — ti) /N. 



6. Summary and outlook 

Motivated by the competing dynamics of different social groups, we analyzed the birth, 
growth and relaxation of a network formed by two populations. The main distinction 
between these two is their preference for having different degrees: Ka and Kb. With a 
society of extroverts and introverts in mind, we considered mainly systems with Kq = 100 
and Kb = 10. In our study, we focus only on the evolution of the network, i.e., the 
connections between individuals. In particular, during each time-step of our simulation 
(defined by the random selection of A^ individuals), we count the degree (k) of each node 
to determine the probability to create or destroy a link (c (fc), 1 — c (k)). In this work, 
we chose c{k) = 1/[1 + [k/ny], which is effectively a Fermi-Dirac function, starting 
at c(0) = 1, dropping monotonically with k through 1/2 over a range controlled by a. 
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Figure 6. Average degree (fc;,) as a function of t for different system sizes N = 
1000 (circles), 2000 (squares), 3000 (diamonds), 4000 (up triangles) and 5000 (right 
triangles). The data are extracted from simulations realized on 500 realizations, for 
A = 0, Ka = 500, Kb = 10 and a = 100. Inset shows data coUapse of (fcf,) and (cb) 
when plotted against {t — ti) /N. 



In this manner, if the degree of the selected node is too low, it creates a link with a 
randomly picked new "friend." On the other hand, if its degree is too high, a random 
chosen link is cut. 

Performing computer simulations, we evolve an initially unconnected set of iV = 
Na + Nb individuals, with a variety of cr's and A^'s. All data collected reveal three time 
regimes, characterized by non-monotonic behavior of the introverts' degrees. In the early 
regime, all individuals are creating links, so that the connectivity in both groups grows 
linearly. The onset of the intermediate regime is signaled by the introverts reaching their 
preferred degree while the extroverts remain unsatisfied. Not surprisingly, the cross-links 
continue to be created by the latter. The introverts struggle to cut links, but, doing so at 
random (in our model), find themselves with more cross-links than a sustainable value. 
In other words, the cross- links overshoot, while the introvert-introvert links undershoot. 
Finally, the late regime is characterized by the extroverts having reached their preferred 
degree, so that the final, stationary state is reached through the redistribution between 
intra-group links and cross-links. To describe these scenarios with more than words, we 
turn to various theoretical approaches. 

A full description of the system can be formulated by a master equation for 
P({/ijj},t), where /ijj = 0//ijj = 1 corresponds to the absence/presence of a link 
between individuals i and j. Though we can write such an equation, we cannot solve it. 
Instead, we consider an approximate formulation, using only three "global" variables: 
a (/3), the total number of links between the extroverts (introverts) and X; the total 
number of cross-links. Even in this reduced phase space, we cannot solve the master 
equation exactly. To make progress, we exploit a standard mean field approach, writing 
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equations of evolution for the average degrees within each population ((fcaa) and (fcbb)) 
and those associated with cross- links: (kab) ■ The emergence of the three time regimes 
can now be predicted and qualitatively good agreement between theory and simulations 
are found. Surprisingly, in the case Na = Nh, the ratio 9o of the restricted degrees 
(kab) and (kbb) is given by an integer in the asymptotic regime of both transient and 
stationary states. Of course, there are limitations to this crude theory, some of which 
we duly noted. 

Our study should be regarded as the first step towards a fully interacting adaptive 
network, in which the individuals are characterized by some states (e.g., opinions, 
susceptible/infected, etc.) as well. In this work, only the links are endowed with 
dynamics. Even so, there are many unanswered questions which deserves being pursued. 
Exploring ways to go beyond mean field theory, by studying the full distribution in k, 
is just an example. Work on this aspect is in progress and will be presented elsewhere 
|48] . Having more than two groups is another. In reality, every individual will have its 
own preferred degree (kj) or flexibility (o"i), which generates an interesting problem with 
quenched disorder. Needless to say, we can study further variations, each incorporating 
more aspects of real social networks. 

If we endow individuals with states and their dynamics, then the vistas are even 
wider. Since one of our motivations for the present study is how epidemics can 
be affected by social behavior, let us provide a simple example by considering the 
epidemiological SIS model on such a dynamic network. See |19] for a review of the 
basis of epidemiological theory and network theory. In this framework, models that 
involve adaptive networks have been considered in [501 EH El], where the networks 
"rewires" by having nodes shifting their connections from one friend to another. If we 
admit the simplest psychology and let the preferred degrees k depend appropriately on 
the level of the epidemic (through the total fraction of infected individuals, say), we 
can expect oscillations to occur. Promising preliminary results have been obtained [53] 
for a much simpler version. In that model, mean field theory predicts only fixed points 
with no limit cycles (though underdamped oscillations can occur). Simulations which 
include stochastics show only a simple stationary distribution, although a signature of 
the feedback induced oscillations can be found in the autocorrelations (and its power 
spectrum [5^155]). In this regard, we believe that novel and interesting phenomena will 
emerge if feedback is introduced into the system by coupling of the dynamics of the nodes 
with that of the links. If realized, such a system will have its own natural frequencies 
(periods) that depend only on the system parameters. How these "internal" periods 
interact with those from "external" drives (e.g., seasonal variations of infectious diseases) 
will surely be of interest in the modeling of epidemics in interacting populations. 
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For our analysis, we consider the A^ — )■ oo limit only, i.e., equations (1231) . (IMl) . (1251) and 
(126|) . For convenience of notation, let us define 

x,y,z= {kaa),{kab),{hb) and /i± = l/(l±A), 
so that 



x + fi+y 



Tb 



z + fi_y 



(57) 
(58) 



Therefore, for arbitrary f (r^) and h (r^), we have 

dxf = f'/Ka, dyf = fi+d^f, (59) 

dzh = h' /Kh, dyh = fi-d^h. (60) 

A useful guide is to keep in mind our specific case, in which — /' = cr/ (1 — /) and 
—h' = ah{l — h) showing explicitly the O (a) nature of these derivatives. Let us write 
the equations for dtx and dtz as 

dtx = Xf — Xi and dtz = Zh — Zi, 



where X = 1 + A + Xi and Z = 1- A + Zi with 



Ai 

For dty, one obtains 
dtV = 



2x 



x + fi+y 



-, and Zi 



2z 



z + i2.y 



(61) 
(62) 



^(4-X)/ + l^(4-Z)/. 



1 + A 



Xi 



1-A 



Zi-2. 



(63) 



Defining the vector 

/ {kaa) - (kaa)* \ 
0(t)= {kat)-{kab)* , (64) 

\ (hb) - {hb)* J 

and linearizing equations ( 1M|) . (l25l) and (l26l) around the fixed point, we write 

ai0(t) = -M0(t) . (65) 

We see that every element of M will contain /' or h' (or both), which are O (a). Denoting 
this order by Mq we have 

/ u Ua-U \ 



V 



u yU+n + fi-V V 
u_f V 



(66) 



where 



U = X*{-f'Y/Ka, U 

V = Z* {—h'Y //tfe, V = 



_ 1 + A 

2 

1-A 



(4 -X*) (-/')* K, 
(4 _ ^*) (_/,y /^,, 
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denote the above functions at the fixed point. In the large a hmit of interest here, 
the simple expectation is that all eigenvalues will be proportional to a. However, this 
would lead to an instantaneous relaxation into the steady state (i.e., with O (1/cr) — )■ 
decay time). A more careful examination reveals that Mq has a zero eigenvalue, with 
eigenvectors 

\ 

(67) 



-/^H 



with 



(0| - ( ^ (1 - 1^) , 1, 2^ (1 - ^) ) and |0) = 1 

(0|0) = |;^ + |^. (68) 

Since both its trace and its co-factor (respectively, v + u + fi+u + fi^v and uv + fi+uv + 
fi^vu) are non-zero, the other two eigenvalues are O (cr). Expecting both to be positive 
on physically grounds, we do not pursue their precise values. 

Given that we have the eigenvectors at O (cr), we can compute the eigenvalue at O (1) by 
standard perturbation theory. Calling that part of the matrix Mi, the 0(1) eigenvalue 
is given by standard perturbation theory: 

A={0|M,|0)^i|:^, (69) 

where Mi is the O (1) terms in M. The decay of any quantity that has overlap with |0) 
will be independent of a in the large a limit. 

Acknowledgements 

We thank S. Dorosz, B. Schmittmann, and L.B. Shaw for illuminating discussion. This 
research is supported in part by the US National Science Foundation through DMR- 
0705152. 



References 



[1 
[2 

[3: 
[4; 

[5 
[6 
[7 
[8 
[9 

[lo; 
[11 

[12 



Steven H Strogatz, 2001 Nature 410 268 

Rcka A and Barabasi A L, 2002 Rev. Mod. Phy. 74, 47 

Dorogovtsev S and Mendes J F F, 2002 Adv. Phys. 51 10791187 

Newman M E J, 2003 SIAM Rev. 45 167256 

Boccaletti S et al, 2006 Phys. Rep. 424 175308 

Costa L d F ct al, 2007 Adv. Phys. 56 167242 

Hubermann B et al, 1998 Science 280 95 

Kleinberg J M, 2000 Nature 406 845 

Adamic L A, 2001 Phys. Rev. E 64 046135 

Bilke S and Perterson C, 2001 Phys. Rev. E 64 036106 

Burda Z ct al, 2001 Phys. Rev. E 64 046118 

Noh J D and Rieger H, Phys. Rev. Lett. 92, 118701 



Network Evolution Induced by the Dynamical Rules of Two Populations 21 

[13] Wah Wu C, 2007 Synchronization in complex networks of nonlinear dynamical systems, World 
Scientific Publishing Co. Pte. Ltd., London 



[14 

[15 
[16 

[17: 
[is: 

[19 
[20 

[21 
[22 
[23: 
[24 
[25 

[26: 

[2?: 
[28: 

[29 
[30 
[31 

[32 

[33: 

[34: 

[35 

[36: 

[3?: 

[38 
[39 
[40 
[41 
[42 
[43: 
[44: 
[45 

[46: 
[4?: 

[48 
[49 
[50 

[51 
[52 
[53: 
[54: 
[55 



Arenas A ct al, 2008 Physics Reports 469 93-153 

Coleman J ct al, 1957 Sociometry 20 253 

Valcnte T, 1995 Network Models of the Diffusion of Innovations, Hampton, Cresskill, NJ 

Kauffman S A, 1993 The Origins of Order, Oxford University, New York 

Keeling M J, 1999 Proc. R. Soc. London, Ser. B 266 859 

Solomonoff R and Rapoport A, 1951 Bull. Math. Biophys 13 107 

Rapoport A, 1957 BuU. Math. Biophys. 19 257 

Weigt M and Hartmann A K, 2001 Phys. Rev. Lett. 861658 

Moukarzel C, 1999 Phys. Rev. E 60 R6263 

Newman M E J and Watts D J, 1999 Phys. Lett. A 263 341 

Newman M E J and Watts D J, 1999 Phys. Rev. E 60 7332 

Moore C and Newman M E J, 2000 Phys. Rev. E 61 5678 

Moore C and Newman M E J, 2000 Phys. Rev. E 62 7059 

Newman M E J et al, 2000 Phys. Rev. 84 3201 

Kuperman M and Abramson G, 2001 Phys. Rev. Lett. 86 2909 

Johansen A and Sornette D, 2000 Physica A 276 338 

Tadic B, 2001 Eur. Phys. J. B 23 221228 

Blatt M et al, 1996 Phys. Rev. Lett. 76 3251 

Zhou H, 2003 Phys. Rev. E, 67 061901 

Reichardt J and Bornholdt S, 2004 Phys. Rev. Lett. 93, 218701 

Newman M E J, 2004 Phys. Rev. E 69 066133 

Newman M E J, 2004 Eur. Phys. J. B 38 321 

Palla G et al, 2005 Nature 435 814 

Derenyi I et al, 2005 Phys. Rev. Lett. 94 160202 

Ravasz E et al, 2002 Science, 297 1551 

Vicsek T, 2002 Nature 418 131 

Guimera R et al, 2003 Phys. Rev. E 68 065103 

Watts D J et al, 2002 Science 296 1302 

Pastor-Satorras R and Vespignani A, 2001 Phys. Rev. Lett. 86 3200 

Pastor-Satorras R and Vespignani A, 2001 Phys. Rev. E 63 066117 

Halloran M E et al., 2008 PNAS 105(12) 4639 

Noh J D et al, 2005 Phys. Rev. E 71 036131 

Boguna M et al, 2004 Phys. Rev. E 70 056122 

, Amaral L A N, et al., 2000 Proc. Natl. Acad. Sci. U.S.A. 97 11 149 

Liu W J, Schmittmann B, and Zia R K P, to be pubhshed 2010 

Keehng M J and Eames K T D, 2005 J. R. Soc. Interface 2 295-307 

Gross T et al, 2006 Phys. Rev. Lett. 96 208701 

Shaw L B and Schwartz I B, 2008 Phys. Rev. E 77 066101 

Schwartz I B and Shaw LB,. 2010 Physics 3 17 

Platini T, to be published 

Alonso D et al, 2007 J. R. Soc. Interface 4 (14) 575-582 

Simes M et al, 2008 J R Soc Interface 5 (22) 555-66 



